X 



DAMTP-HEP-92-24 

NSF-ITP-92-75 

May 1992 



Mark Hindmarsh 



^ ■ Semilocal Topological Defects 

Department of Applied Maths and Theoretical Physics 
t^^ I University of Cambridge 

Cambridge CBS 9EW 
U.K. 

> . 
0\ , and 

(N 

(Nl ' Institute for Theoretical Physics 

University of California 
CN ■ Santa Barbara CA 93106 

0\ : U.S.A. 

ix. 



Abstract 

Semilocal defects are those formed in field theories with spontaneously broken 
symmetries, where the vacuum manifold M is fibred by the action of the gauge 
group in a non-trivial way. Studied in this paper is the simplest such class of 
theories, in which M ~ S"^^^^, fibred by the action of a local U{1) symmetry. 
Despite M having trivial homotopy groups up to 'K2N-2^ this theory exhibits a 
fascinating variety of defects: vortices, or semilocal strings; monopoles (on which 
the strings terminate); and (when N = 2) textures, which may be stabilised by 
their associated magnetic field to produce a skyrmion. 



1 Introduction 

The spontaneous breaking of symmetry lies at the heart of modern particle physics, for 
if we believe that nature is more symmetric at higher energies, then we must also have 
an explanation for the apparent lack of symmetry in the universe today. The success of 
the partial unification of the electromagnetic and weak forces in the framework of the 
Salam- Weinberg gauge theory |^ has led to many proposals for enlarging the symmetry 
groups of matter and forces, perhaps most famously Georgi-Glashow SU(5) unification 
0. Global symmetries have also been suggested, notably the axial and family |Q 
ones. When we combine these ideas with the paradigm of the hot Big Bang |^, finite 
temperature field theory tells us that the universe would have gone through a series of 
phase transitions at which more of the symmetries would have become hidden or 
broken (although there are proposals for symmetry increase with reduced temperature 
0). In the past decade or so it has become clear that the effects of the phase transitions 
can be felt long after they finish, through the formation of topological defects R M 



and massless Nambu-Goldstone modes [T^. In fact, these phenomena are almost the 
only way of deriving information about the unification scale of 10^^ GeV. 

Topological defects, which may result from the spontaneous breaking of global or 
local symmetries, can briefly be characterised as regions of space-time with dimen- 
sion d < 4 in which the order parameter of the phase transition vanishes due to the 
topological properties of the broken and unbroken symmetry groups. In decreasing 
order of their dimension, they are termed domain walls [^, strings or vortices |^, ^, 
monopoles [ll[], and finally textures [|12|, |T^ or sphalerons [jl4[. These last are single 



spacetime points where the order parameter goes through zero. Related to this last 
class are skyrmions [^ , which can be thought of as textures prevented from collaps- 



ing by higher derivative terms in the action. Nambu-Goldstone modes result from the 
breaking of global symmetries only, and they seem to possess an important property 



which allows them to persist long after the phase transition is complete: scaling ||T0|. 
Because they are massless, the only relevant length scale for their dynamics in the early 
universe is the Hubble length H~^. The direction of the order parameter is correlated 
over this distance, so at any epoch the energy density in the Nambu-Goldstone modes 
is ~ rfH'^, where rj = |0|. The density parameter for these modes is then roughly 
Grf^ where G is the gravitational constant. This scaling property is also thought 
to be shared by strings (both local and global) ||T^, global monopoles [|1^, and global 
textures [|13]. If this is the case, then any of these objects will produce a scale-free spec- 
trum of density perturbations. It is tempting to regard it as more than a coincidence 
that the required amplitude for the formation of structure in the universe is obtained 
if 7] is about the Grand Unification scale. The energy in the Nambu-Goldstone modes 
and in topological defects may also account for the recently observed fluctuations in 
the Cosmic Microwave Background [|T8]. Other observable effects include gravitational 
lensing (principally from strings) and the production of a significant density of radia- 
tion, either gravitational or Nambu-Goldstone. The radiation density is bounded by 



the expansion rate of the universe during nucleosynthesis |jT9[, in the same way that 



the number of neutrino species is bounded (see e.g. Kolb and Turner 0). 

The above considerations therefore motivate the study of topological defects, par- 
ticularly strings, monopoles, and textures, and also global symmetry breaking, in a 
cosmological context. A piece of ideology picks out the class of models considered 
in this paper: it is that nature should accomplish its phase transitions in the most 
economical way, with the minimum of scalar fields. The scalar sector of most GUTs 
so far proposed is made complicated enough in the effort to reduce the gauge symme- 
tries. This leads one naturally to the study of models in which both global and local 
symmetries are broken at the same phase transition. 

In this paper, the simplest class of models which exhibit this property is studied. 
One might term them extended Abelian Higgs (EAH) models, for they consist of N 
scalar fields with a U{N) symmetry, but only with the overall phase gauged, leaving 
an SU{N) global symmetry. When A^ = 1 we have the ordinary Abelian Higgs model. 
It has already been shown that this class has some unusual properties, and it contains 
some general lessons for the study of phase transitions and the associated topological 
defects. Perhaps the most remarkable point is that the models exhibit stable vortices 
for all A^, despite the first homotopy group of the vacuum manifold being trivial when 



A^ > 1 |2^, |21[. This is made all the more interesting by the observation made by 
Vachaspati [^ that the N = 2 case represents the scalar sector of the Electroweak 
theory in the limit that the Weinberg angle is 7r/2, and that the vortices retain their 
stability away from this limit. The stability of the vortices of the EAH models was 



examined in Ref [^, and it was shown there that they are stable only if a certain 
parameter /5, the ratio between the scalar self-coupling and the square of its gauge 
charge, is less than or equal to 1. The case /5 = 1 is particularly interesting, for here 
the equations of motion can be reduced to a first order form, the Bogomol'nyi equations 
23| . In this case it was shown in Ref |21[] that instead of there being a unique vortex 



solution (the Nielsen-Olesen vortex EM for A^ = 1), there is a family of them labelled 



by a complex (A^ — l)-vector. Furthermore, it was shown in Ref. |^ that there exists 
a whole family of n-vortex solutions labelled by Nn complex numbers, the extra 2n 
real parameters being identified with the positions of the centres of the vortices. 

In Section 2 of this paper, the model is introduced and the stability of the vortices 
demonstrated by numerical methods. It was pointed out in [^ that these vortices 
bear a very close relation to the 2-dimensional (2D) lumps in CP^^^ a-models p6 |. 
In Section 3 the connection is made in more detail, it being shown that the long- 
wavelength dynamics of the EAH models are those of the cr-model. When /3 > 1, the 
next term in an expansion in field gradients is shown to be a Skyrme term [|T3] . 



In 3 space dimensions there are further unusual features to the EAH models. It was 



surmised in pO[ that because the vacuum manifold is simply connected, the vortices 



could end in what was described as a "cloud" of field gradient energy. In |^ it was 



asserted that this cloud was in fact a sort of global monopole, but for reasons of space 
the field configuration was not written down. In Section 4 this defect is remedied, 
and it will be shown that the structure resembles a Dirac monopole with a physical 
string. The energy of the configuration is of course infinite, as is that of a single global 



monopole |T^ or a infinite length of string. 



Tfiere is also another 3D solution of interest in the special case N = 2 where the 



vacuum manifold is a 3-sphere: the texture |12, 13|. This is explored in Section 5, 
where it will be seen that it possesses an associated magnetic field in a remarkable 
configuration best described as a twisted vortex. In the /? ^ oo limit, there is the 
possibility that the magnetic field pressure can prevent the final collapse and unwinding 
of the scalar field. The magnetic field produces a fourth derivative term in the effective 
(T-model lagrangian - the Skyrme term in fact - and thus the resulting object can be 
called a skyrmion []T3|. 



The question of the formation of all the above objects at a cosmological phase tran- 
sition, and their subsequent evolution, is tackled in Section 6. The average magnetic 
flux through a correlated area can be computed exactly in the N = 2 case, and from 
this it is shown that, even when they are stable, vortices are rather rare. They are also 
likely to be short, terminating at monopoles. However, when /? < 1 the attraction of 
monopoles to antimonopoles is argued to be greater than the string tension, so there is 
a mechanism for the creation of string through the annihilation of a pair of oppositely 
charged magnetic poles attached to different segments of string. It is then plausible 
that a cosmological scaling solution can be set up, in which at any epoch the Hubble 
volume contains a few collapsing textures, annihilating pairs of poles, and segments 
of string. When j3 > 1 the dynamics resemble more that of a a-model, since it is 
favourable for the scalar field to remain in its vacuum manifold. In either case, a scale 
free spectrum of density perturbations would be established, although their statistics 
would probably be quite different. It has already be proposed that, separately, string 
P, ^, global monopoles ^^, and global texture |jl3| could seed the structure that we 



see today. The model here presented can produce all three types of objects, at rela- 
tively little cost in complexity, so the notion that it may be of importance to cosmology 
is not far-fetched. 

In order to back up this notion, some justification for thinking that such models 
have some basis in a realistic particle physics theory. To this end, a Grand Unified 
model possessing the required characteristics is presented in the concluding Section 7, 
along with some discussion of the results obtained in this paper. 

2 Extended Abelian Higgs models and their vortex 
solutions. 

The model studied in the bulk of this paper is a theory of iV complex scalars $ with their 
overall phase gauged and an SU{N) global symmetry. The most general renormalisable 
lagrangian in 4D consistent with these symmetries is 

C = \D^<^\' - IXm' - r^Y - \f,,F^'' (1) 

where D^ = d — ieA^ and |$p = $$. At zero temperature the field has a vacuum 
expectation value of magnitude 77, and the symmetry of the vacuum is reduced to a 



global U{N — 1). In particular, the gauge symmetry is broken, and the physical spec- 
trum consists of a vector particle with mass rriy = ^/2er], 2{N — 1) Nambu-Goldstone 
bosons, and a Higgs scalar with mass rus = v2A?7. The symmetries of the potential are 
stable to radiative corrections, because |$p is the only invariant consistent with the 
SU{N) X f/(l) invariance. The issue raised in Refs. [pO| , |2T1| was whether there are any 
non-trivial classical solutions to the equations of motion following from (|ID. It is clear 
that there are when A^ = 1, for then the theory reduces to the ordinary Abelian Higgs 
model, for which it is known that there are relativistic flux tubes or Nielsen-Olesen 
vortices |Q. In fact, any classical solution of the Abelian Higgs model carries over to 



the EAH models. This can immediately be seen from the equations of motion 

D2$ + A(|<l>p-r/2)$ = 



(2) 



which are solved by any fields of the form (A^, (/)$i/?7), where $i is a constant unit 
vector, if (A^, 0) are solutions of the ordinary Abelian Higgs model. However, it 
does not necessarily follow that they are stable solutions, for it may be that a small 
perturbation away from the embedded Abelian Higgs model solution grows. From 
another point of view, one might not expect there to exist vortex solutions at all, for 
the conventional indicator of the presence of vortices, the first homotopy group of the 
vacuum manifold M, is trivial for all A^ > 1. However, we shall see that this fact is to 
be interpreted with some caution, because (as shown in |^) stable vortex solutions, 
termed "semilocal strings" in Ref. pO], do exist, but only when f3 = A/e^ < 1. 



The simplest vortex solution is of course a static straight string, which can be 
chosen to lie on the z axis with fields independent of z. Accordingly, this dimension 
will be ignored for the moment, and the transverse dimensions shall be denoted x* with 
i G {1,2}. Searching for static solutions is equivalent to minimising the energy (per 
unit length) functional 



S= (fx 



\DM^ + \\m^-V^? +\B^] (3) 

where B = eijdiAj is the {z component of) magnetic flux. If S is to be finite, each term 
in (2.3) must seperately vanish faster than r^ at infinity. The vanishing of the potential 
energy density implies firstly that |$| — > 77, i.e. that $ lies in M at r = 00. The most 
important observation from our point of view is that the covariant derivative term must 
also vanish at large r, which implies that $(a;*) must be a gauge transformation: that 
is, 

Hx') -. exp (ze£ dy'Aiy)] $(a;^) (4) 

when |xo| and |a;| tend to infinity. The continuity of $ demands that the integral of 
the gauge field around the z axis be an integer multiple of 27r/e. This quantisation 
condition divides the finite energy field configurations into topologically inequivalent 
sectors, but says nothing about the existence of vortex solutions, nor does it guarantee 



that they are stable. The magnetic flux must also disappear as r -^ oo, which imposes 
a boundary condition on Ai, that it must become pure gauge sufficiently fast. 

We will ffist examine solutions with cylindrical symmetry, for which we may choose 
a gauge in which 

$ ^ $ie^"i^ Ai -^ meijXj/er^ (5) 

For simplicity we shall restrict ourselves to solutions in the rii = 1 sector, for which 
the most general solution has the form 



N 



* = E UO^ce'""-' ci^ = e^,^AO/e (6) 

Q = l 

where ^* = er]x'^ are dimensionless coordinates, giving distances in units of ^/2 times the 
vector particle's Compton wavelength. For the rest of this section we shall also rescale 
the scalar field $ -^ ^rj so that it becomes dimensionless. The boundary conditions 
on the fields are 

fa -^ Sal, a — > 1 as ^ ^ oo 
a — > as C "*■ 

and the fa vanish at ^ = if tIq, 7^ 0, which includes of course /i. The imposition of 
cylindrical symmetry means that the effect of any spatial rotation around the z axis 
can be undone by a global U{N) transformation. This forces the $„ to be orthogonal, 
and we may write 

<l„<l>^ = r]^5af3 (7) 

As Vachaspati and Achucarro showed [^, there exist solutions with all the fa except 



/i vanishing, which are essentially Nielsen-Olesen vortices with winding number ni. 
We denote the solutions /i and a in the special case rii = 1 by / and a. Near the 
origin / ~ ^ and a ~ .^^, while the behaviour at infinity is 

/ ^ 1 - cir^ exp{-^0 « ^ C2e^ exp(-v^O (8) 

with Ci,C2 constants P3. Thus both the scalar and vector fields tend exponentially 



to their vacuum values outside regions which are of order m~^ and m~^ in width 
respectively. 

The question raised earlier was whether these solutions are stable or not. In order to 
answer it we must expand the energy functional to quadratic order around the solution 
to find the fluctuation operator, and then look for negative eigenvalues. Only if there 
are no negative eigenvalues is the solution stable. Generally, this operator is not diago- 
nal in the perturbations a^, (p, as there can be terms of the form / d'^C,{S'^S/6Ai6^)aiip 
in the expansion. However, we can remove these terms by a suitable choice of gauge, 
one developed by Moss, Toms and Wright p7|: 



ditti - i^i^ip - ip<^) + Uitti = (9) 



where the real vector f/j satisfies the equation 

t/,|<l>|2 = 15,|$|2 (10) 

In this gauge, the quadratic part of the variation in the energy functional becomes 

lA^r + pm' - 1)1^1' + iPiv'^ + '^v?] (11) 

where D is the covariant derivative with respect to the background field Ai. When the 
background field is an embedding of the Nielsen- Olesen vortex, $^ = (/e*^, 0, . . . , 0), 
^i = ^ijijO'li'^1 the quadratic variation splits into two pieces: one in Oj and the upper 
component of $; and the other in the remaining components, or 

5^£ = 5'£no + E (lAv^ar + PiP - ^)\Va\') (12) 

a>l 

The subscript NO denotes that this part of the variation is precisely that obtained from 
the ordinary Nielsen-Olesen string, which we know to be stable against all perturbations 
tti and (fii. Therefore we need consider perturbations in the lower components of the 
scalar field only. Because of the global symmetry we may chose the perturbation in any 
direction with no loss of generality. Thus writing ipa = Va J2n V^nC*""^ with Va constant 
vector, the question of stability reduces to solving the eigenvalue equations 



■i|(4)+^(-«)^+«f-i; 



^n = eipn (13) 



If the perturbation operator in square brackets has any negative eigenvalues, then there 
will be growing modes, which means that the Nielsen-Olesen string is merely a saddle 
point of the energy functional S. 



Firstly, we establish that there is at most one unstable mode. Equation (|l^) has 
the form of a Schrodinger equation with potential f/„ = (a — n)^/C,^ + f3{f^ — 1). Now, 
if Un > Um for all ^, then the same inequality applies to the lowest eigenvalues €Q{n) 
and eoim). Therefore if the lowest eigenvalue with potential [/„ is positive, it follows 
that for those m satisfying Um > Un for all ^ there are no negative eigenvalues. We 
know from the stability of the Nielsen-Olesen vortex that Ui produces no negative 
eigenvalues, so stable perturbations will be obtained if 

(a-n)2> (a-l)2 Va G [0, 1) (14) 

Thus |t2| > 1 give stable perturbations, leaving one possible unstable mode, which has 
n = 0. 

I have been able to make further analytic progress only for (3 = 1, which shall be 
described later. For the rest, therefore, I resorted to numerical methods. To solve 
(O) one first must find the background fields / and a. To do this I used a simple 



relaxation method [^. The energy functional is discretised onto a ID lattice (taken 
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to be 800 points), with derivatives evaluated on the hnks, and then minimised at each 
point r in turn. This means solving a cubic equation for fr and a linear one for a^ to 
obtain their updated values. One can also incorporate an overshoot factor to speed 
up the convergence. For example, if /„ is the solution to dS/dfn = 0, then /„ is 
updated to f'n = fn + c^ifn — fn), with a > 1. The energy converges reasonably 
fast: typically, after a few hundred iterations it is changing by a factor of less than 
10"^. The accuracy of the algorithm can be checked by calculating E/2'Krf for the 
Nielsen-Olesen vortex with (3 = 1, for it should be exactly 1 (as will be explained 
below). Using FORTRAN double precision with a step length A^ of 1.25 x 10~^ the 
program gives a value 0.999993. The resulting / and a were then fed into (|13D which 
was then solved numerically by a shooting method. The differential equation (11) 
was discretised by a simple symmetric scheme whereby ip'^ = {ipr+i — '^r-i)/2A^ and 
ip" = (ipr+i + i^r-i — '^^'r) / ^^'^ ■ Negative eigenvalues are found by shooting for ipj. = at 
the boundary r = Vmax, with initial conditions ipo = 1 and ipi = 1. A Newton- Raphson 
iteration on the function ^iprmax (^) "W'as used to automate the convergence. One can find 
the lowest eigenvalue by a judicious choice of the first two values of e, and checking 
that the eigenfunction has no nodes. The results for /? = 1, 3, 10, 30, and 100 are 
listed in Table 1. For P = 0.3, 0.1, 0.03, and 0.01 no negative eigenvalues were found. 
The accuracy of the eigenvalues is good: they have been quoted to 3 decimal places, 
because doubling the step length and halving the number of points changed the lowest 
eigenvalue for /5 = 10 from -9.1654 to -9.1653. The accuracy is not quite as good as 
this suggests, because as we shall see next the lowest eigenvalue for /5 = 1 is not about 
4 X 10^'^ but is exactly zero. 

In order to prove the existence of this zero mode we need a bit of machinery, 
established by Bogomol'nyi. He showed that in the special case /3 = 1 the second order 
field equations split up into first order ones. This emerges from rewriting the energy 
functional ^: 

+ |(/5-l)(l$P-l)^] (15) 

where ni has been taken to be positive, and all currents are assumed to vanish at 
infinity. At /3 = 1 it can be seen that the energy is minimised at 2Tmiiif when the 
following first order equations are satisfied: 

{Di + it^,D,<^)^ = (16) 

5+|$|2-l = (17) 

In the cylindrically symmetric case with ni = 1 the Bogomol'nyi equations are 

fL + lia-n^)fa = a' + e(E/'-l) = (18) 

From these equations it is easy to check that the zero mode of the vortex is 



ijj = A/'exp 



f 



(19) 



where A/" is a normalisation factor. This zero mode is a symptom of an unexpected 
degeneracy in the solutions to the Bogomol'nyi equations. 

Without loss of generality we may take the zero mode to be in the $2 direction. 
Note that the 3 Bogomol'nyi equations are not independent, for we may replace /2 by 
fiw/^ with any complex w and still satisfy the equations of motion, providing /i and 
a solve 

/{ + ^^y^/i = «' + ^[/i'(i + kl V^') - 1] = (20) 

li w = we are left with the Bogomornyi equations for the plain Abelian Higgs model, 
for which we know there exists a vortex solution. In Ref. [^] it was proved that vortex 



solutions to ([T7| ) exist for any complex w, revealing the Nielsen-Olesen vortex as just 
one of a family of vortex solutions. Recall now that the choice of $2 is arbitrary so 
long as it is orthogonal to $1 and has unit modulus. This reveals that the family is 
labelled by a complex (A^ — l)-vector w = w^2- The asymptotics of this family for 
w y^ are very different from those of the Nielsen-Olesen string, whose fields approach 
their vacuum values exponentially (H). Expanding /i, /2 and a up to fourth degree in 
( = ^/w it is easily found that 

/i ^ i-lr^ + d-H-'K-' 

/2 ^ C'-lC', (21) 

a ~ 1 - C^ + {I - A\w\~^)C^ 

while the magnetic field B tends to 2\w\^^(^'^ = 2|u7p^^^. This power law behaviour 
is quite strange: the width of the flux tube is an arbitrary parameter instead of the 
Compton wavelength of the vector boson. This is not at all what we expect for magnetic 
fields in the Higgs phase of a theory. 

It is clear from (|2l| ) that |w| controls the width of the vortex. The magnitude of 
the scalar field at the origin also depends on this quantity, and it is possible to derive 



the bounds 21 



l-2|w|"^ < /2(0) < 1 (22) 

showing that as the width of the vortex tends to infinity, the magnitude of $ gets closer 
to 1 everywhere, even at the vortex core. 

The existence of general solutions to the Bogomol'nyi equations (0) has been 



proved in Ref. |2^, where it was shown that in the n- vortex sector they are labelled by 
nN complex numbers. These numbers can be interpreted as (A^ — 1) internal complex 
parameters for each vortex, plus another n labelling their positions. There are also time 
dependent solutions, corresponding to charged vortices, with their own Bogomol'nyi 
equations [^. 

Clearly, the Bogomol'nyi limit f3 = 1 is rather special, so we are entitled to ask 
what happens to the vortices away from this value. Let us consider an isolated vortex 
whose field configuration satisfies the Bogomol'nyi equations. Its energy functional can 
be rewritten 

£ = 2nr]^ + \{I3 - l)ri' j (f^B^ (23) 

8 



using the second of equations (p!?!). Recall that at large distances, B — > 2\w\'^/^'^: hence 
the second term in (|23|) can be expressed as a function of w, or i(/5 — l)r7^fi:(|w|)|w|~^. 
For large \w\, n tends to a constant, because the dependence on \w\ is merely an 
expression of how the magnetic field energy changes under a scale transformation. 
Thus, according to the sign of (/3 — 1) the energy will change with the size of the 
vortex. The tendency is to collapse if /3 < 1 and expand if /3 > 1 , which is consistent 
with the perturbative stability at w = 0. In the strict Bogomol'nyi limit it is possible to 
approximate the low energy dynamics of the field theory by slow motion on the space of 
parameters, or moduli, of the non-trivial classical solutions ||5^ . For the single vortex, 
omitting the trivial centre-of-mass motion, the moduli are the real and imaginary parts 
of w. The time dependence of the fields occurs entirely through the moduli changing 
with time, so that the infinite dimensional lagrangian for the fields 

de[|$|V + |(^.)']-^ (24) 

reduces to a finite dimensional one in w. In the one vortex sector it is logarithmically 
divergent, since $ contains a piece w/i/^. We can cut off this divergence at a scale R, 
which might be interpreted as a vortex seperation, to find 



v' 



T{w)\w\^ - U{(3 - l)\w\-^] (25) 



where T{w) includes at piece proportional to ln{R/\w\). From this we can compute the 
timescale for the collapse or expansion of the vortex outside the Bogomol'nyi limit. If 
we just consider the radial motion in the moduli space, and ignore the time dependence 
of the logarithm, the solutions to the equations of motion are 

\w\^ = \w{to)\^ + ait - hf\w{h)\-^ (26) 

where a = k{j3 — l)/2ln{R/\w\). Thus the timescale for changing size by a factor 2 is 
approximately w(to)^/\/«- We have been working in units where times and distances 
are measured in units of m~^, so if W is the width in natural units, the collapse 
timescale goes as W'^niy. Thus vortices which are large compared with the gauge 
boson Compton wavelength collapse very slowly, in much more than the light crossing 
time. 



3 The CP^ ^ cr-model as the low momentum hmit 

In this section it is shown how the low momentum dynamics of the EAH models with 
A^ scalars approach those of the CP^^^ a-model [^, where by low momentum we 
mean momenta whose magnitudes are much less than the masses m^ and rriy. Let us 
first describe the finishing point: cr-models with target space CP"^. These manifolds 
are the set of all complex lines in C™^^: thus Z'^ and Za in (7™-+^ are the same point 
in CP"^ if and only if there exists a complex number a such that Z'^ = aZ^- When we 



use homogeneous coordinates Za on (7"^+^ the natural induced metric on the imbedded 
Qpm jg ^]-^g Fubini-Study metric [plf 



9ab{Z) = ^^ (27) 

This metric is used to construct the a-inodel action from m + 1 complex scalar fields 

ScP^ = J d'x (d^^agabmO^'^b) (28) 

It is easily found that this action is invariant under multiplication of the fields by a 
space-time dependent complex number, and so it has both a local scale invariance in 
the target space and a local f/(l) gauge invariance. To arrive at equation ( pS] ) we first 
rewrite the gauge potential 

The expression for the gauge field arising from this field redefinition is 



F =-- 

^•^ 2e 



a(^$9,)<l>|$|^-(9(^$<D)($a,)$) /\W + f,u (30) 



where f^u = di^^a^). Complex geometers will recognise the first term on the right hand 
side as 1/e multiplied by the components of the Kahler form uj^^ in the immersion $ : 
i?" -^ (^m+i 1^ After some algebra the lagrangian of the EAH model is reexpressed 
as 

- Y^^,.r - \f,.r - \Km' - vr m 

The degrees of freedom are now explicitly separated. There is a real scalar field |$|, the 
Higgs mode, which because of the potential term tries to keep as close to rj in magnitude 
as possible. This gives a mass \/2eT] to the propagating modes of the gauge field a^, 
which in turn react back on |$|, with the effect of reducing its size. The first term in 



(31) can be immediately identified as the (7P™ a-model model lagrangian, describing 
the NG bosons, provided $ remains close to 77. Again, spacelike excitations in the $ 
field tend to reduce |$|, but provided the typical momentum p in these excitations has 
IpP -C Xt]^ — ml, they will have little effect. The back reaction from the gauge field is 
slightly more involved. Consider the equation of motion for a^: 

a^r'^ + 2e2|$|2a" = --df.uj^" (32) 

The right hand side acts as a source for a'^, and is third order in derivatives in $. Thus 
in a momentum expansion we can say a^ = T]0{p^/{er]Y). The back reaction of the 
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gauge field on the Higgs mode can therefore be neglected if {erjY ijr' / {erff'Y <^ Xrf, or 
p^/ml < 13. 

A term not appearing in the usual a-model action is that arising from the gauge 
field, —u^yOj^^/Ae^ = 0{p^/e^). It is consistent to retain this as the next highest order 
term if it is larger than (|$p — ?7^)(9^$a(yfab5^$fe). The equation of motion for the Higgs 
mode tells us that deviations of |$p from t]"^ are of order p^/X. Therefore the condition 
to be met is 0{p'^/e'^) <^ 0{p'^/X), or /? -C 1. A major result from this section has 
now been arrived at: the low momentum effective lagrangian for the EAH model in 
the (T-model limit (/3 <^ 1) is 



1 
4?' 



£eff = V^d^^aQabd^'^b - T^^,..^^' (33) 



Recall that the Kahler form u;^,^ is the first term on the right hand side of (^) multiplied 
by 2e. There is a special significance to the additional term in (|33|): it is effectively 
a Skyrme term |]r3| for the CP"^ cr-model. This is because it is fourth order in all 



derivatives but only second order in derivatives with respect to time. There is yet 
further geometrical significance to the relationship between the gauge field and the 
Kahler form u. Firstly, uo^ /m\ is the volume form on CP"^, where by u"^ we mean the 
wedge product taken m times [0. For maps $ : i?^™ -^ CP'"^ there is a topological 
winding number Q^-, whose expression in local coordinates is 

Q„= 1 / ,-,,„i|%|%...|t|iL) ,,.......,.„ (34) 

This winding number is an integer for fields tending to a constant at spatial infinity 
in i?^™, and it measures the number of times $(x) covers the target CP"^. These 
integrals are just the Chern classes 0]. Q2 is the 4D abelian anomaly, while Qi is the 
total flux in the plane, in units of 27r/e. Thus for the CP^ cr-model, the flux through 
any 2-surface S measures the area of the target space covered by the image of $. In 
higher dimensional target spaces the relation between the restriction of the Kahler 
form to $(S) and the induced volume form dV breaks down. This is a consequence 
of Wirtinger's Inequality, which states that for a 2k-dimensional submanifold Sk of 
M, uj^/k\ < dVk, [0. The equality holds only if $(5'^) is a complex submanifold of 
the target space, which is a severe restriction. For 2D submanifolds it amounts to 
demanding that $„ can be written as f{x){l,wi, . . . ,Wm), with the Wa holomorphic 
(or antiholomorphic) functions of a complex coordinate z = x + iy on H. 

Our effective lagrangian (^) reveals that although the manifold of minima M has 
the topology of 5*^^"^, as far as the global degrees of freedom are concerned it is 
really CP^~^. This is the result of the existence of a special set of circles in M, those 
generated by gauge transformations. Every point in M lies on such a circle, and in 
fact this set of circles forms a fibration of S*^^"^, the famous Hopf fibration |^. The 
canonical example occurs at A^ = 2, for CP^ ~ S*^, and it is well known that S^ can 
be constructed as a twisted bundle with fibre S^ and base space S'^ (or vice versa, of 
course). The physical significance of these fibres is that the scalar field can lie along 
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them without costing any gradient energy, for any derivatives can be counterbalanced 
by a gauge potential This is the underlying reason for the existence of semilocal 
strings, for the condition of finite energy per unit length forces the field to lie on a 
Hopf fibre at infinity. The number of times n this fibre is traversed is given by 



1_ [ dx^^^ = — 

inJsL l$P 27r 7^1 



n = -- L dx^^ = - L dx^A, (35) 



The second equality follows from the fact that Oj must vanish at infinity for finite energy 
solutions. The last term is of course nothing but the total magnetic flux through the 
plane in units of 27r/e, which we saw previously was equal to the number of times the 
image of $ overs the target CP^. We now see explicitly that these winding numbers 
are one and the same, even away from the low momentum limit. The integral of 
the magnetic flux is \ j d?x{e~^ujij + fij)€ij. For finite energy solutions the second 
term must vanish because of the boundary condition imposed on Oj at infinity. This 
leaves an expression which may be integrated by parts to reveal 2TTn/e. One may 
legitimately object that Uij is not defined at points where $ = 0. However, it is 
defined for every neighbouring field configuration in which $ vanishes nowhere, and 
for these configurations the equality holds between the fibre winding number and the 
degree of the map from R^ to CP"^. Thus we can consistently extend the definition of 
Qi to maps $ which vanish somewhere. 

One of the principal reasons for the interest in CP"^ cr-models is that they possess 
instanton solutions in 2D, which have an arbitrary size due to the conformal invariance 
of the theory. The analogy with 4D Yang-Mills theories, coupled with the nice prop- 
erties of 2 dimensions, make these cr-models ideal laboratories for instanton physics. 
In fact, as we shall see, these instantons are just semilocal strings in the large scale 
limit. The CP"^ instanton solutions can be written down very simply in closed form. 
A single instanton is ^E\ 

« = '^-^°'*'+'"'*' (36) 

[\z - za\^ + \w\f 

where Zq and w are arbitrary complex parameters and $i$2 = 0. This solution rep- 
resents an instanton centred at Zq with size \w\. The similarity to the EAH model 
vortices in the Bogomol'nyi limit is especially clear if we expand the instanton field at 
large distances and compare with (pl|): 

$ = $1— 1 - 1^ + -^ + $2^ 1 - \-^ 37 

\z\ \ ^ \zY 4 \z\^ j \z\ \ ^ \zY j 

Note that at infinity the field tends to $ie*^, which is a single point in CP^, and 
identical to the vortex configuration. The differences only start to appear at 4th degree 
in |w|/|^| in this expansion. Perhaps the biggest difference is that vortices exist in the 
EAH model even at w = 0, where the CP"^ lump is singular. This singularity causes 
problems in numerical simulations of the scattering of CP"^ lumps |^B[ , and so the low 
energy scattering of semilocal vortices in 2D should be far better behaved |^ . 
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To conclude this section we specialise to the case N = 2, where the vacuum manifold 
is S^ and the a- model target space is CP^ c^ S"^. This isomorphism is realised through 
the Pauli matrices a^, which enables us to construct three real fields 

^A ^ 1,^^$/^ (38) 

The effective lagrangian ( |33D can be rewritten in terms on these scalar fields, and is 
1 



This may be verified using the relations c^ta'^f^ — '^^ab^ha ~^ ^aa^bb^ ^^^ imposing the 
constraint |$p = f]"^. Equation (|39| ) is the lagrangian for the 0(3) cr-model with 
a Skyrme term. We shall see in the next section that the equivalence between the 
(T-models is only local and can break down in an interesting way in more than 2 
dimensions. 



4 Monopoles in Extended Abelian Higgs Models 

The fact that tt2{CP"^) = Z naturally leads to the question of whether monopoles 
exist in the theory. Since the coordinates on CP"^ are ungauged degrees of freedom, 
these would have to be a type of global monopole. The properties of the EAH model 
monopoles are best illustrated in the case m = 1 through the equivalence to the 0(3) 
(j-model ||2^, which is well known to possess singular global monopole solutions |]ll 



If, however, we view the constraint (p^cp^ = rf as being imposed by a potential, i.e. 

Coiz) = Id.^^'d'^^^ - ^A(0V - VY (40) 

then the theory possesses well-behaved monopole solutions of the form 

0^ = nx^hir) (41) 

where r^ = x*a:*, /i(0) = and h{oo) = 1. The energy of these solutions is linearly 
divergent, because the energy density at large distances is simply ?7^/2r^, due to the 
angular variation of the field. (In the gauge monopole []TT| the angular derivatives are 



cancelled by the gauge potential.) Thus if i? is a large radius cut-off, the energy of a 
monopole is 

Em = 2nRr]^ + E^ (42) 

where E^ is the contribution from the centre of the monopole at which |0| —>■ 0. The 
energy divergence makes the discussion of a single global monopole slightly problematic, 
and dependent on boundary conditions. For example, if the centre of the monopole is 
kept fixed, there is a cylindrically symmetric zero mode which takes the monopole away 
from the hedgehog configuration of ( |1T| ) p5[ . There is in fact a 1-parameter family of 



azimuthally symmetric monopole field configurations, all with the same energy (^ 
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A physical cut-off would be imposed, for example, by a neighbouring antimonopole 
at a distance ~ R. Since dEM/dR is independent of R, the force between the pair of 
poles can be argued to be generally independent of their separation (although it also 
depends on the details of the field between them ||3^), which makes them very efficient 
at annihilating in the early universe W7 . 



What would such a configuration look like in the CP^ coordinates $? It is easy to 
check that the field configuration which gives (ETl) at large distances is 



/sini^e*^\ ,^^. 

where {6, (j)} are polar coordinates. However, this configuration is ill-defined at ^ = tt. 
In CP^ there is no sign of this problem, since (e*"^, 0) is just a single point. There are 
two ways of resolving the issue. The first is to treat it as a gauge artifact, and define 
$ in two different coordinate patches covering the upper and lower hemispheres of S"^ 



oo 



(the sphere at r = oo) [|ri|]. Equation ( ^31) would be the field on the upper hemisphere 
$(7, while the field on the lower would be 

$ . = ( ^7 i^ . ^ (44) 

It is easy to see that ^u and $i are related by a (singular) gauge transformation where 
they overlap. The other resolution of the problem is to look for solutions in which $ 
vanishes at = tt in order to avoid the singularity. With this departure from spherical 
symmetry the full solution should have the form 

\ h2{r,U) cos ^9 I ^ ' 

where /ii = = /i2 at r = and hi{r,7r) = 0. As we approach the z axis near 
6 = n and at large r the fields should be more or less z independent. Thus the top 
component of $ should behave more or less like the field of a Nielsen- Olesen vortex, as 
it obeys similar field equations and changes phase by 2n around the axis. The bottom 
component vanishes as p'^/r'^, where p^ = r"^ — z"^. Thus we find a string on the —z axis. 
To support this interpretation, let us compute the gauge potential and the magnetic 
field. At large distances, away from the string, $ should lie in its vacuum manifold 
and its covariant derivative should vanish, which means that 

sin^i^ . ^ 

A, = 0, ^ 46 

er sm t/ 

where 0i is a unit vector in the azimuthal direction. Near 6 = n, A^ blows up as 
l/(ep). This is precisely the field that would be produced by a string on the —z axis, 
with fiux 27r/e in the upward direction. The magnetic field resulting from (^) is 
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This is a rather strange result. We started by writing down a global monopole config- 
uration: from it the EAH model has produced a local monopole, with one Dirac unit 



of magnetic charge 27r/e [11|. This remarkable object is in fact the 6w = 7r/2 limit 
of a field configuration first singled out by Nambu in the Electroweak theory, which 
is essentially a string carrying Z^ fiux terminating on a (electro)magnetic monopole 
Il37| . He proposed that an approximately stable classical solution might be obtained 
by spinning a finite length of string with opposite poles at each end. 

If the string at —z is viewed as a gauge artifact, in other words as a Dirac string, 
then we have constructed a truly spherically symmetric field configuration, albeit at the 
price of violating one of Maxwell's equation, diBi = 0. Consistency can be maintained 
either by invoking General Relativity to hide the source of fiux in a black hole, as 



has been investigated by Gibbons et al. [^, or by embedding the theory in one with 
a larger non-Abelian symmetry, so that the source of the fiux is a 't Hooft-Polyakov 
monopole. 

In considering the string/monopole system we are again faced with the problem 
of infinite energy common to all global monopoles. The energy function, cut off at a 
distance R, can be split into three components: the monopole core, the string, and the 
long range monopole field. The energy can then be written 

E = E^ + 2r(Tfu{(3)R + 2'KprfR + -n/e^R^ (48) 



where z^(/3) is a monotonically increasing function, having the value 1 at /3 = 1 
The third term is the energy in the scalar field gradients, with p a number depending on 
the exact field configuration, but which is known to be bounded below by 1 |^| . This 



bound is saturated when all the scalar field gradients are concentrated on a line starting 
at the pole. The fourth term is the energy in the magnetic field, with R^. the monopole 
core radius. For sufficiently large R the magnetic field energy is negligible. We can 
try to circumvent the infinite energy problem by considering a monopole-antimonopole 
pair, whose strings run along the z axis from —R/2 to — cx) and from R/2 to -|-cx) 
respectively. Although this configuration has infinite energy also, it is meaningful to 
talk of the difference in energy AE between it and a straight string on the axis. Naively, 
this energy difference is just twice the energy of a monopole with cut-off i?, minus the 
energy of the missing piece of string, or 

AE ~ 2Ec + 2TipTfR - 2txu{I3)R (49) 

Thus for Type I vortices {(3 < 1) the energy is minimised when R = 0, whereas for Type 
II vortices (/3 > 1) the monopole separation can increase without bound, provided a 
suitable configuration is chosen for the global degrees of freedom. This behaviour is 
consistent with the results concerning the stability of the vortices as a function of (3 
in Section 2. The creation of a pair of poles in a straight string is essentially a very 
large perturbation in the scalar field away from the Nielsen- Olesen string in the region 
—R/2 < z < R/2. We saw that at /3 < 1 extended non-singular vortex configurations 
tended to collapse towards the Nielsen- Olesen solution, and it is encouraging that 
this heuristic analysis indicates in another way the reappearance of the Nielsen- Olesen 
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string. Numerical work p9[ also seems to confirm the tendency of non-singular vortices 



to collapse or expand according to the sign of /3 — 1. 

5 Textures and Skyrmions 

To complete the study of topological defects in EAH models, we consider two related 
types of defects that occur when the third homotopy group of the vacuum manifold, n^, 
is non-trivial. They are textures |TB[ and skyrmions [|T3| , which are field configurations 
which map the whole of 3 dimensional space, with the point at infinity included, onto 
a non-contractible 3-sphere in the vacuum manifold M. A simple example is the 
0(4) Goldstone model, whose Lagrangian is essentially (|iO| ) but with 4 real fields. A 
spherically symmetric field configuration which covers M once is 

01 = cos X <p2 = sin X cos 6 

03 = sin X sin 6 cos ip 04 = sin x sin 6 cos Lp 

where xl'"? t) vanishes at r = and reaches vr at r = cxd. Exact solutions to the a-model 
equations of motion were found by Spergel and Turok [^ , which collapse from infinite 



size at t = — oo to a singular point and then rexpand towards t = +oo. At the point 
of collapse, the gradient energy density becomes sufficiently large to force the scalar 
field through = 0, which removes the topological winding. The dynamic nature of 
this object is a necessary consequence of Derrick's scaling argument in 3D |^^, which 
shows that both the gradient and potential terms in pO|) decrease with scale, and so 
any localised object can decrease its energy by shrinking. If, however, we supplement 
the energy functional by a term which is fourth order in derivatives. 



gM<Pn\d,rY + g2{d,<P^d,<P^Y\ (51) 

then the possibility arises of a stable minimum in the energy for configurations such 
as (^), because |ii^4| increases linearly with decreasing scale. If gi+ g2 = 0, the action 
remains second order in time derivatives, and E^^ is positive. This favoured choice is 



the so-called Skryme term |T^ . The importance of skyrmions lies in their providing a 
low energy model of the nucleon, where the Goldstone bosons are interpreted as pions 
and the topological winding number as baryon number. 

The vacuum manifold of an EAH model with N scalars is 5^^^^, and so the only 
one with textures or skyrmions has N = 2. (It is well known that 7Tk{S"') = for 
k < n [^.) We have already seen that the low momentum dynamics of the model are 
equivalent to the CP^ a-model, which is in turn equivalent to the 0(3) a- model. Thus 
we can expect the large textures of the N = 2 EAH model to behave in a very similar 
way to the textures in nematic liquid crystals, since their effective theory is a S"^ / Z2 



cr-model [^. In particular, we cannot expect spherical symmetry. 

In the cr-model approximation to the theory, we can parameterise the constrained 
scalar field by angles C, "0, and a as follows: 

/ siniCe^('^+^)/M ..o^ 
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It is easy to show from 0"^ = ^a^^/rf that C, and if) are polar coordinates on the global 

CP^ target space. The gauge potential that results from ensuring that the current 

vanishes everywhere is 

1 



A„ 



2e 



[cos C^^li' 



^^.cr) 



showing explicitly that a is just a gauge mode. The gauge field is then 



flU 



2e 



sinC(9(^^9^)C) 



(53) 



(54) 



s ^ e 



At this point it is instructive to recall the remarks made in section 3 concerning the 
geometrical interpretation of F^^ as being (up to a factor of e) the volume form on the 
image of $ in the target CP^~^. In 2D, it was shown that the flux though a surface 
S was proportional to the area of the $(S). Here, we can display the equality very 
clearly 

/ miC,dC,dil} (55) 

/*(S) 

The right hand side is easily seen to be the area on the target space for the global 
degrees of freedom, S"^. 

In studying textures and skyrmions in this theory, it is not entirely clear which 
coordinate system to use, for one can no longer assume spherical symmetry for the 
objects. Skyrmions in the 0(3) a-model have been studied by de Vega and were 
rediscovered later by Wu and Zee WM- They both exhibited toroidally symmetric 
configurations which covered the target S"^ exactly once, namely 

C = C(/i) i^ = (3-^ (56) 

where {/i, /?,0} are toroidal coordinates ^^. The boundary conditions are C(0) = 
and C(oo) = vr. However, this cannot be the basis for a toroidally symmetric solution, 
because in these coordinates the second order part of the energy functional is 



E2 = 2'KRrf / dnd(5 



sinh/i 



(cosh n — cos /9) 



'dC 

dfi 



+ sin^ ( coth^ jj 



(57) 



where R is the radius of the circle /i = 00. The explicit dependence on f3 shows the 
failure of the toroidal ansatz, which includes the spherically symmetric case through 
the stereographic projection cosC = 1 — 2sech^/i. 

Nonetheless, it is illustrative to compute, in the low momentum limit, the magnetic 
field that this configuration produces. Defining r to be cosh fi — cos (3, we find 



A, 



B, 



2eR 

T 



2ei?2 



cosC(A 



^i 



T d ( sinh n 



sinh fxdn \ 



T 



COS C + 0iT 



dfi 



T 



■cosC 



(58) 
(59) 



(Some useful formulae to aid this computation can be found in Ref. [0.) The lines of 
fiux encircle both the z axis and the circle fi = 00 |^^. The fiux crossing any surface of 
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constant (3, which is a spherical bowl with /i = oo as its border, is 27r/e, while the flux 
crossing surfaces of constant y?, which are half planes with the z axis as boundary, is 
also 27r/e. Thus this field configuration can be interpreted as a twisted vortex ring. In 
fact, the lines of flux trace the Hopf fibres in the projection S^ -^ R^ that this texture 
embodies. 

There are two topological invariants which characterise the textures and Skyrniions. 
There is the Hopf number P^ , ^, ^, which measures the degree of the map $ : 



i?^ U oo ~ S''^ — > 5^, and is given by 

^^ = A I d^'xe'^^^l^.^d^^agabdk^b/m^ (60) 

Stt^ J 

where Qah is the Fubini-Study metric on CP^~^. There is also the Chern-Simons 
number 



e2 



^-=2^ 



Jd^x^'^AF^k (61) 



In the low momentum limit, where the gauge field is determined entirely by the scalar 
field, these two quantities are equal, as equations (PD| ) and (^) quickly show. For the 
EAH model texture Nh = Ncs = 1- 

The toroidal ansatz is particularly good for picking out the way the energy behaves 
under a change of scale. The energy in the magnetic field is 



TV 

E4 = I d^d(3- 



Ae'^R J sinh /x 



d f sinh /i 

d^\ T 



cosC 



2 



+ sinh /i 



'cosC' 



2N 



(62) 



E4 clearly behaves in the opposite way to E2, showing that the Skyrme term can 
stabilise the texture in the cr-model approximation, (3 -^ 00. Comparison of the relative 
sizes of the second and fourth order terms shows that the skyrmion has size of order 
l/erj, implying a mass of order r]/e. Away from the limit, it becomes possible for the 
scalar field to be forced through $ = if its gradient energy density gets too large. It 
is therefore an open question whether a skyrmion exists in the N = 2 EAH model for 
finite /3. 

A quantum skyrmion will at best be merely metastable, since the field can tunnel 
through the barrier which prevents unwinding, namely the potential V{^). The Eu- 
clidean action for such an event is of order V{0){eri)~'^, and so the skyrmion decay rate 
is estimated to be 

r, ~ {eri)e-^''/'' (63) 

where 7 is a numerical coefficient. In view of the possibility of factors of 27r appearing 
from the symmetries of the skyrmion, 7 need not be merely of order unity. 

If there is no stable skyrmion then the unwinding results from the classical evolution 
of the collapsing texture [^. In either case, the magnetic field will then dissipate, as 
there are no longer any scalar field gradients to support it, and the Chern-Simons 
number changes from 1 to 0. If there are any anomalous fermionic currents coupled 
to our model, there will be a violation of charge conservation associated with the 
unwinding event. This process can be identified with the manufacture and decay of a 
sphaleron in the Electroweak theory in the limit 6^ = tt/2. 
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6 Cosmological Phase Transitions 

We now turn to the possibility of the creation of any or all of the defects described in 
the preceding sections. It was pointed out by Kibble [^ that when the cooling universe 
goes through a phase transition the order parameter (which here is $) is uncorrelated 
beyond some length scale ^. At a second order transition, this scale is determined by the 
correlation length of the scalar field at the temperature at which thermal fluctuations 
cease to explore the vicinity of $ = 0, or the Ginzburg temperature [^. At a first 



order transition it is the average separation of the bubbles which nucleate the transition 
which sets the scale. Thus different domains in space are mapped to different regions 
of the manifold of minima of the finite temperature effective potential, Mt- If this 
manifold has non-trivial homotopy groups ttq, tti or 1^2 then it is possible to produce 
defects (walls, strings and monopoles respectively) at the boundaries of the correlated 
domains, where the continuity of the order parameter forces it to vanish. 

One can estimate the probability of defect formation by triangulating both space 



and Mt, and then assigning points in Mt at random to points on the spatial lattice ^_ 
This simulates the appearance of domains of size ^ at the phase transition. When M^ 
has simple geometry, such as that of a sphere, it is possible to calculate the probabilities 
without triangulation, just by assigning points on the sphere to the vertices of the 
spatial lattice and interpolating geodesically between them [Q. This method seems 



more accurate as it does not involve a drastic truncation of the field configuration 
space. For example, let us recall the calculation of the probability per unit correlation 
volume of forming monopoles, where the target manifold is S"^. Firstly, we triangulate 
Mt with a tetrahedron. A monopole is formed if the boundary of a 3-simplex in the 
spatial triangulation, which is of course a tetrahedron, is mapped onto the tetrahedron 
approximating Mt- There is then a topological obstruction to mapping the interior of 
the spatial 3-simplex into Mt-, which means that the order parameter must leave Mt 
and vanish somewhere inside. This is the monopole. The total number of different 
ways of mapping the vertices of the spatial triangulation into Mt is 4^ = 256. The 
number of these mappings which are mappings of one tetrahedron onto another is just 
24, the order of the tetrahedral group, including reflections for antimonopoles. By this 
method, therefore, the probability of forming a monopole or an antimonopole is 24/256 
= 3/32. In the alternative method [^], the image of the vertices of a spatial 3-simplex 
forms 4 random points on the target manifold. The gradient energy density of the 
field is minimised if the vertices are joined by geodesies, and the edges filled in by the 
smaller of the two possible spherical triangles. A monopole is formed if the sphere 
is covered by this construction. To evaluate the probability of covering the sphere, 
consider any three of the four image points, plus their connecting geodesies, extended 
to great circles. These circles divide the sphere into 8 spherical triangles. The sphere 
is covered by the boundary of the spatial 3-simplex if (and only if) the fourth point 
is contained within the triangle antipodal to the one defined by the original 3 points. 
The probability of this happening is proportional to the average fractional area of a 
spherical triangle. By this construction the probability is seen to be 1/8. 
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In the case at hand the situation is rather more coniphcated. We wish to assess the 
formation probabihty of three types of object: strings, monopoles, and in the N = 2 
case, texture. For the first two, the relevant topological quantity is the magnetic flux 
through a surface E. If the surface is closed, the Maxwell equation diBi = tells us 
that the total flux through the surface is zero. There are no true magnetic monopoles 
in this theory, only approximately spherically symmetric objects which are the termini 
of strings. If the surface has a boundary, the total flux JF tells us how many strings 
(more strictly, the difference between the numbers of strings and antistrings) could pass 
through that surface. This is just the integer part of eJF/27r, say n. However, unlike 
the ordinary Abelian Higgs model, there is no topological condition which forces the 
scalar field to vanish somewhere on this surface, since its manifold of minima is simply 
connected. Nevertheless, when (3 < 1, the configuration near the surface S with least 
energy will consist of a set of Nielsen- Olesen type vortices with total winding number 
n, plus any stray flux left over. In a statistical ensemble like that produced at a phase 
transition, however, the most likely configuration is the one with the minimum free 
energy, so it is not necessarily correct to conclude that the flux is everywhere confined 
in vortices. In fact, there is a lot of entropy in short lengths of string terminating on 
monopoles, so it seems probable that the phase transition produces string segments 
rather than infinite string. 

Let us now compute the flux through a surface E with dimensions of the correlation 
length ^. As usual, we triangulate space with a simplicial complex whose average 
edge length is ^, and choose S to be one of the 2-simplices of the complex. In using 
the topological properties of M^, we are assuming that the gauge field is entirely 
determined by the gradients of the scalar field, which becomes an increasingly good 
approximation in the limit e — >■ cx). In this limit thermal fluctuations cannot generate 
real gauge bosons as they are too masssive. We now assign values of the scalar field, 
constrained to lie in Mp, at each vertex. The field is extended to all points in R^ by 
joining the images of the vertices in Mt with the shortest possible geodesic. Since M^ 
is simply connected there is no topological obstruction to then filling in the resulting 
set of geodesic line segments with geodesic surfaces, so that $ is defined everywhere 
on S. The task of computing the magnetic flux through this surface is simple when 
N = 2, for as we saw in Section 3 it is proportional to the area of the image of the 
simplex, as measured by the Fubini-Study metric on CP^. Now CP^ is isomorphic 
to 5*^, so the average fractional area occupied by the image of a spatial 2-simplex is 
1/8, exactly as for the monopole calculation. If the surface were to cover the whole 
of the target CP^, as in the vortex, the total flux would be 27r/e. Hence the average 
flux through a 2-simplex is 7r/4e. The flux through neighbouring 2-simplices is by 
construction uncorrelated, so that means we need an average of 8^ = 64 contiguous 
2-simplices before collecting enough magnetic flux to form a vortex. Then when f3 < 1 
and vortices are stable, we expect the string density to be suppressed by a factor of 8'^ 
relative to the A^ = 1 case |^. 

When there are more than 2 components to the scalar field the calculation is not so 
straightforward, because the integral of the magnetic field over a 2-simplex is no longer 
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equal to the area of its image in CP^^^. As was discussed in Section 3, by Wirtinger's 
inequality the flux is almost always less than the area. 

Let us now turn to the probability of texture formation. First we reprise the cal- 
culation for pure global texture on a vacuum manifold isomorphic to S^. The simplest 
triangulation of Mt is into the boundary of a 4-simplex, a 4D tetrahedron, which has 
5 3-simplices, 10 2-simplices, 10 edges and 5 vertices. What is required is the average 
topological charge per correlated volume: in other words, the average volume in the tar- 
get S'^ per spatial 3-simplex. This gives the number of texture collapses that will occur 
as the universe expands, for it appears from dynamical simulations |^ that a region 
with charge greater than 1/2 [i.e. covering more than half the 3-sphere) will eventually 
collapse and unwind. There can only be charge on a spatial simplex if it maps onto 
a 3-simplex in M^- The probabilty of this is 5!/5^ = 24/125. The 5! assignments are 
divided into two equal sets with topological charge ±1/5, so the average charge is of 
course zero, but he average of the magnitude of the charge is 24/625 ~ 0.0384. The 
alternative calculation |4^ gives the average volume (and thus the average topological 



charge) of the 3-simplex defined by 4 random points as 2~^ = 0.0625. 

For texture in the N = 2 EAH model, we have to exercise care in assigning meaning 
to the topological charge, because of the gauge invariance in the model. The gauge 
invariant quantity which measures the volume of the 3-sphere is the Hopf number, 
defined in equation (^). One can check |^ that the integrand, which we term the 



Hopf density H, really is the volume element on 5''^ by writing $-^ = (yi +iy2, ys + iyAJ. 
Imposing the constraint |$p = 1 it is found that 

2 

T-i = —eijkdiyidjy2dky3 (64) 

1/4 

The average topological charge density is now computed, as before, by assigning random 
values in S^ to the vertices of a spatial 3-simplex. As before, it is 1/16. 

The evolution of an EAH model after a phase transition seems to be a rather 
complex problem, especially for N = 2. For large f3 the scalar field is forced to stay in 
Mt and we need not worry about the formation of strings and monopoles as they are 
unstable. The dynamics are those of a CP^~^ a-model, which should be tractable in 
the large-A^ limit. If 0{N) a-models are anything to go by [|l^, it seems reasonable 
to suppose that a scaling solution should be established, with the correlation length 
of $ increasing in proportion to the Hubble length H~^. Accompanying the scalar 
field there will be a long range magnetic field belonging to a broken the broken gauge 
symmetry, a peculiar prediction of this class of models. However, the strength of this 
field is negligible, of order H'^/e'^, and it contributes only a tiny fraction ~ G^H^/e^ to 
the energy density of the universe. There is an interesting complication when N = 2, 
for here there are textures, a few of which should collapse per Hubble volume per 



expansion time [13|. If /? is sufficiently large, a population of skyrmions may result, as 
the magnetic field pressure prevents the final unwinding of the "knot" of scalar field 
gradients (see Section 5). Skyrmions, even at the GUT scale, may represent a serious 
cosmological problem, for their lifetime is ~ ]^Q-39g7/3/47ra g (^^ ^ ^^d we are perforce 
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discussing the large j3 case. For example, if 7/3 > 0(10) then the lifetime becomes 
longer than Is, which brings with it the potential for conflict with standard Big Bang 
nucleosynthesis 0. 

At low /?, the stability of the vortex solutions opens up the possibility that they 
dominate the evolution of the model. Let us assume, as argued above, that the phase 
transition results in the formation of finite length segmenta of string terminated by the 
monopoles of Section 4. If these were local monopoles one would expect the strings 



to disappear in a very short time [11] as the string tension pulled the monopoles 
together. However, there is a competing effect: the linear potential between monopole- 
antimonopole pairs, which can act to join two nearby segments into a single longer 
one. This string can be created by the annihilation of global monopoles. It is thought 
that both global monopoles and strings, separately, can reach a scaling solution in an 
expanding universe, so it is conceivable that the (3 < 1 EAH model can also scale. 
At any time, each Hubble volume would contain a few lengths of string, terminating 
on annihilating global monopoles. The string would decay in the usual way by loop 
production and gravitational radiation. 

7 Discussion and Conclusions 

We have seen that EAH models exhibit a surprisingly rich set of classical solutions. 
They also teach some valuable lessons about defects in field theories. When there are 
spontaneously broken local and global symmetries, the standard test for defects, the 
existence of non-trivial homotopy groups of the vacuum manifold, must be applied with 
care. Let the group of local symmetries of the theory be Gi, which must be a subgroup 
of the group of all the global symmetries G, since any local symmetry is also a global 
one. The action of the gauge group on the vacuum manifold divides it into a set of 
orbits, each isomorphic to the coset space Gi/Hi, where Hi is the unbroken subgroup 
of Gi. Let us call this coset space Mi. The action of choosing a gauge picks out a 
single point on each orbit, and the set of these points forms a manifold Mg, which is 
the target manifold of the Nambu-Goldstone fields. What we have, of course, are all 
the elements of a fibre bundle [^. The bundle space is M, the fibre is Mi, and when 
we choose a gauge we project down to the base space Mg. If it were assumed that 
M is globally a direct product Mi x Mg, the homotopy groups of M would then be a 
direct sum 7ik{Mi) + iik^Mg), and any non-trivial topology in Mi and Mg would show 
up in 7rfe(M). As we have found, however, this is not a safe assumption to make in 
general. If we had made it for the N = 2 EAH model, we would have concluded that 
M c::^ S^ X S"^, and that the non-trivial homotopoy groups were a local tti, and a global 
712 and TTa. This would indicate that this theory had gauge vortices, global monopoles, 
and global texture. The reality is somewhat more complicated: the gauge vortices are 
not necessarily stable; the global monopoles have a magnetic charge supplied by a real 
string attached to it; and the global texture also has an associated magnetic field. This 
is all a result of the fact that it is possible to assemble a fibre and a base space to 
make a smooth manifold in more than one way: in this particular case, M is the Hopf 
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bundle [O, with fibre S"^ over a base space S'^, which is globally isomorphic to S^. In 



view of this feature I offer a definition of a semilocal defect: a defect in a theory whose 
vacuum manifold is a non-trivial bundle with fibre Gi/Hi. 

As a result, one must reconsider the statement sometimes made that d-dimensional 
defects exist in n space dimensions if 7r„_d-i(M) is non-trivial. Suppose we are search- 
ing for finite energy defects. The finite energy condition means that, at spatial infinity, 
the field must lie in its vacuum manifold and be on a gauge orbit: that is, it must lie in 
Ml. Thus one must really examine T[n-d-i{Gi/ Hi). However, not even the non-triviality 
of this group guarantees the existence of a stable static solution to the field equations, 
for we saw that the stability of the vortex solution in the EAH models depends on the 
sign oi (3 — 1. Thus a non-trivial homotopy group reveals only that finite energy field 
configurations fall into inequivalent classes, and does not generally give information 
about the existence of stable solutions within these classes. 

At the start of the paper, I put forward arguments in favour of models with mixed 
local and global symmetries, based on a principle of economy in Grand Unified The- 
ories. The class of models considered here have turned out to have much intrinsic 
interest, but it would be encouraging to demonstrate that GUTs really can produce 
semilocal defects. Serious model-bulding is beyond the scope of the present work, a 
candidate semilocal grand unified model is at last possible. Let us try to embed the 
local symmetry in a plausible grand unified group, which must therefore be of at least 
rank 5. The most obvious group to try is 50(10), and its simply connected covering, 
S'pm(lO). The Nielsen-Olesen string can be embedded in this theory by first break- 
ing it with an adjoint to SU{b) x f/(l). This f/(l) may now be broken with a 126 
dimensional 0i26 down to Z2, resulting in stable strings [^]. It is of course possible to 
break directly to SU{b) x Z2 (strictly, 5f/(5) x ^10/^5 0) with the 126 alone, but 
the resulting string is not simply an embedding of the Nielsen-Olesen vortex [^. The 
126 contains an SU{5) singlet $ which accomplishes the required symmetry breaking, 
and it is this component that makes the string. In order to give a global symmetry to 
the string field we must therefore replicate the $126- 

Finally, the reader will have noticed that no attempt has been made to calculate 
the cosmological perturbations produced by these models. This is because the calcu- 
lation seems rather messy (except perhaps in the limit of large N and large (3 where 
the (T-model approximation works well). There exist calculations in various stages of 
refinement for strings pBI and texture (Spergel and Turok |T3) separately, but it is 



not clear how they are combined with the global monopole signal in EAH models. It 
is however worth pointing out that existing texture calculations of CMB fluctuations 
assume a vacuum manifold isomorphic to S''^. A collapsing texture in an A^ = 2 EAH 
model will not give the same distinctive signal, a 10° hot or cold spot on the microwave 
sky, if only because of the lack of spherical symmetry. This may be of importance to 



the texture galaxy formation scenario JTSl, which may be in danger of conflicting with 



the recent COBE data flS 
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8 Tables 

TABLE 1: Table of eigenvalues e of equation (|I3]), for various values of the parameter 
(3 (see text). 



/5 


e 


100 


-99.234 


30 


-29.208 


10 


-9.165 


3 


-2.088 


1 


3.81 X 10-3 
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